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Abstract 

A scenario of galaxy formation is put forward which is a process of sud- 
den condensation just after recombination. It is essentially based on the 
fact that the cosmic matter gas after recombination is a general relativistic 
Boltzmann gas which runs within a few 10 6 years into a state very close to 
collision-dominated equilibrium. The mass spectrum of axially symmetric 
condensation "drops" extends from the lower limit M ~ 1O 5 M to the upper 
limit M ~ 1O 12 M . The lower limit masses are spheres whereas the upper 
limit masses are extremely thin pancakes. These pancakes contract within 
a time of about 2, 5 • 10 9 y to fastly rotating spiral galaxies with ordinary 
proportions. In this final state they have a redshift z ~ 3. At an earlier 
time during their contraction they are higly active and are observed with a 
redshift z ~ 5. 



PACS numbers: 98.62.Ai; 04.20. - q 
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1. Introduction 

The formation of galaxies and large-scale structures in the universe is one of the unsolved 
problems in the standard model given by the Friedmann-Robertson- Walker-mo dels. Al- 
though these models are altogether very successful and therefore widely accepted no way has 
been found hitherto to explain in this framework the formation of the observed structures in 
a satisfactory manner. 

The essential reason for this lack is the extreme smothness of the young universe. This 
smoothness has been verified by the COBE data of the cosmic background radiation. Of 
course, COBE observed anisotropics AT/T ~ 10~ 5 of the background radiation only on a 
rather large scale of 7° to 10°. Hence one might observe on a much smaller scale much larger 
anisotropics. Nevertheless this would be rather surprising. So it seems that one has to resign 
to a very smooth universe when it was young. Carrying over the anisotropy AT/T ~ 1CT 5 of 
the background radiation to the density contrast Sg/g ~ 10 -5 at the end of the recombination 
period one is confronted with the problem that these contrasts can not grow up to galaxies by 
gravitational instability in a sufficiently small time [1]. Indeed, today galaxies are observed 
with redshifts z ~ 5 which corresponds to about 10 9 y after the big bang. A further problem 
almost never mentioned is the fact that the galaxies with upper limit masses (M ~ 1O 12 M0) 
are mainly flat and fastly rotating discs whereas the objects at the lower limit (M ~ 1O 5 M ) 
are spheres. 

Today it seems to be proved that the overwhelming part of matter in the universe is non 
luminous dark matter. This is confirmed by the flat rotation curves in galaxies as well as by 
gravitational lensing of galaxy clusters. Many cosmologists believe that this dark matter is 
non baryonic exotic cold dark matter which is subject only to gravitational interaction. They 
propose a scenario of structure formation which is first of all structure formation in the cold 
dark matter, the luminous baryomic matter then adds to these structures. The extremely 
cumbersome numerical calculations yield interesting structures which look very much like to 
those observed in the universe [2]. However, there are at least two problems in this scenario. 
First of all, which particles form the cold dark matter? Hitherto there is no experimental 
evidence for such exotic particles. Secondly, the scenario might perhaps explain the very 
large-scale structures, it seems, however, not to explain the origin of the single objects - 
namely the galaxies and the smaller objects down to the globular clusters. 

In the following paper we propose a scenario of galaxy formation in an universe without cold 
dark matter, that means, the dark matter is purely baryonic. A first version of this scenario 
has been published by the author in two previous paper [3,4]. This version describes only 
the final state of galaxy formation, not, however, the process. This process has now been 
understood, and thus the complete scenario can be represented. 

The basic idea of the scenario is the fact that the cosmic matter gas is from the end of 
the recombination period a general relativistic Boltzmann gas which runs within a few 10 6 
years into a state very close to collision-dominated equilibrium. Hence we sum up briefly 
from the previous papers [3,4] in sections 2 and 3 the properties of the cosmic matter gas 
after the recombination period and the properties of a general relativistic Boltzmann gas in 
collision-dominated equilibrium. In section 4 we describe the process of galaxy formation 
which is a process of sudden condensation. In section 5 we investigate the time evolution of 
the condensation '"drops'". Finally we discuss our results in a conclusion. 

2. The cosmic matter gas after recombination 



-2- 



When the universe has cooled down to T ~ 5 • W^K the recombination of the cosmic plasma 
to neutral H- and He-atoms begins. It ends at 

T ~ 3 • W 3 K = T rek (1) 

when all matter consists of neutral atoms. This recombination period is throughout a drastic 
event in the history of the universe. One of the important consequences is the drastic change 
of the range of interaction between the particles. It decreases from the long range Coulomb 
interaction between charged particles before recombination to the very short range collision 
interaction between neutral atoms after recombination. In an earlier paper [5] the author 
argued that the range r c of collision interaction in the cosmic matter consisting of 70 % 
hydrogen and 30 % helium is approximately given by 

r c ~ 7a (2) 

with Bohr's radius do- Hence we have after recombination the relation 

r c n 1/3 < 3 . W -7 ( 3a ) 

with the particle density 

n = rin + n-He = 0, 8-^- (36) 

(g = mass density, m p = proton mass). From (3a) we can conclude that the cosmic matter 
after recombination is a Boltzmann gas. Futhermore the author calculated in [5] the ratio 
tc/texp where t c is the mean time between two collisions and t exp is the characteristic expansion 
time of the universe. This ratio takes the value 

texp V / 

From this ratio we conclude that the cosmic matter runs from the end of recombination 
(T = T re k) into a state which is very close to equilibrium. But how does this process run 
and how much time takes it? To answer this question we consider in the following section 
first of all the final state of this process, namely the general relativistic Boltzmann gas in 
collision-dominated equilibrium. 

3. The general relativistic Boltzmann gas in collision dominated equilibrium 

We note that the cosmic matter after recombination runs into a state very close to equilibrium 
which is not exact the equilibrium. Clearly, the expansion of the universe does not admit an 
exact equilibrium but only a state very close to equilibrium. Later on we shall assume that 
the physics of this state is essentially that of exact equilibrium. Thus we deal in this section 
with a general relativistic Boltzmann gas in collision-dominated equilibrium. For simplicity 
we assume a one-component gas with particle mass m. 

We adopt the essential properties of a general relativistic Boltzmann gas in collision-dominated 
equilibrium from the comprehensive reviews given by Stewart [6] and Israel [7] . A space-time 
filled by such a gas has a time—like killing field 



(5) 
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where the u l are the components of the 4-velocity and k is Boltzmann's constant (latin indices 
running from 1 to 4, signature +++-). A reference frame /Co can be found such that 

ft = (0,0,0, const) (6) 

In /Co the matter is in rest and the following relations hold 

9ik = g lk {x l ,x 2 ,x 3 ), T = T{x\x 2 ,x 3 ) (7) 

544T 2 = const (8) 

2 

mn(x x , x 2 , x 3 ) = q(x 1 ,x 2 ,x ?i ) = const -T ■ K 2 {--=-) (9) 

rC-L 

where n is the particle density and K2 the modified Bessel function. The specific relation (9) 
between the spatial changes of density and temperature is due to the fact that in a general rel- 
ativistic Boltzmann gas in collision-dominated equilibrium the chemical potential u divided 
by T is spatially constant. We note explicitly that (9) is not an equation of state. It de- 
scribes exclusively the spatial changes of density and temperature in an equilibrium solution 
of Boltzmann's equation. 

In the range of low temperatures kT « mc 2 an asymptotic expansion of K2 yields 

g = const- T 3 / 2 -exp[-^] (10) 

This extremely sensitive exponential relation between the spatial changes of density and 
temperature has been found also by Dehnen and Obregon [8] in a quite different way. 

The metric is determind by the field equations. We consider small anisotropies of the tem- 
perature, that means 

T = T 1 [l + r(x 1 ,x 2 ,a; 3 )] 

(11) 

Ti = const, ^(x^x 2 ,^ 3 )! << 1 
Furthermore we assume the gravitational field to be weak in the sense 



9ik = m + hik{x l ,x 2 ,x s ) 

T) ik = diag(l, 1,1,-1) (12) 

\hikix 1 , x 2 , x 3 )\ « 1 

Starting from eqs. (11) and (12) the author has given a detailed analysis of the linearized 
field equations (with zero cosmological constant) in his previous paper [4]. From this analysis 
follows the basic equation 



At = -kq < ^£ic 2 exp 



mc 2 



fcTi 



(1-r) 



+ aTf (13) 
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where A is the Laplician, kq = 8irG/c 4 Einstein's gravitational constant and a = 7,564 • 
10~ 15 erg -cm^ 3 -K~ A the black-body constant. Eq. (13) is for given background temperature 
T\ and some start value g\ of the density a partial differential equation for the temperature 
anisotropics r in a general relativistic Boltzmann gas in collision-dominated equilibrium. The 
solution r determines then the density by means of 



qi exp 



mc 

~kJ\ 



(1-r) 



(14) 



Eq. (13) says that in a general relativistic Boltzmann gas in collision-dominated equilibrium 
anisotropics of the temperature must exist. Indeed, with r = it follows from eq. (13) that 
f?i = and T\ = and hence g = and T = 0. The existence of these spatial fluctuations of 
temperature and density (by means of eq. (14)) is an intrinsically general relativistic effect. 
It is due to the coupling of the Boltzmann gas with the space-time metric via the field 
equations. In a special relativistic theory, however, space-time is fixed so that field equations 
do not exist. Then we have 544 = const = — 1 and hence T = const and g = const according 
to eqs. (8) and (9). Because of mc 2 jkT\ « 1 (later on we have m = m p and T\ ~ 10 3 i"T so 
that m p c 2 /kTi ~ 10 10 !) it follows from eq. (14) that smallest changes of r lead to extremely 
drastic changes of g. Hence space-time is filled with sharply edged gas clouds although the 
temperature is nearly isotropic. This is the final state in the condensation process described 
later in section 4.2. 

We are interested in a single gas cloud. Such a cloud is described by the equations (see ref. 
M) 



T c Ay = - 



8irG 1 



-g c c exp 



mc 2 T c 



(i-y) 



+ aT 4 



y = -<h y c = y(0,0,0) = 1 



(15a) 



Qc exp 



mc 2 T c 
kTi 



(i-y) 



(156) 



Here the centre of the cloud is placed in the origin of the spatial coordinate system, and g c 
and r c are the values of g and r in the centre. 

Again we emphasize that the density-temperature relation (10) or (15 b) holds in the 
final quasi-equili- brium state . The dynamical process ending in this final state will be 
treated below in section 4.2. 

4. Galaxy formation 

4.1 The cosmological background 

Before describing in section 4.2 the galaxy formation as a sudden condensation process we 
lay out the cosmological model in which this process runs. In fig. 1 the age t u of the universe 
depending on the present day density £0 and the Hubble-parameter Hq = h ■ km/secMpc 
is plotted where k = marks the flat universe. The underlying cosmological model is a 
standard model with zero cosmological constant. 

The recent controversy about the Hubble-parameter between Sandage/Tammann (small 
value of h) and Freedman et al. (large value of h) seems to end in a mean value h ~ 70. 
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Furthermore the earlier large age of the globular clusters of about 16 • W 9 y has been recently 
revised to (12 ± 2) • 10 9 y [9, 10]. So assuming a density £>o,L — 2 • 10 _31 <7 • cm -3 of the lumi- 
nous matter and taking into account the dark matter by a factor 5 to 10 we have a present 
day density go ~ 10 _30 <7 • cm" 3 . From fig. 1 we then read an age of the universe between 
12 • 10 9 y and 13 • 10 9 y which is compatible with the age of the globular clusters. So all is right. 
However, a flat universe as predicted by an inflationary scenario in the very early universe 
seems to be ruled out. Obviously we are living in an open universe. 

The line element of an open standard model in quasi-cartesian coordinates reads 



da 2 



£(x«) 2 



{dx i f + (dx 2 y + (dx 3 ) 2 




(16) 



(greek indices running from 1 to 3) . From the field equations with zero cosmological constant 
we obtain then in the matter-dominated era (T < T re k = 3 • 10 3 K) 



Ki IK 
ct = R (! + -)>---,„ 



(1 + K/R)2 + 1 



(1 + K/R)* 



(17a) 



with 



K = ^hqc 2 qoRoi Rq = ( H o - ^ K oc 2 £o) 2 



(176) 



and where the index o denotes the present day values. In the temperature range T>10 3 -ftT 
equation (17a) can be well approximated by 



ct 



1 



2K2 



From this equation we obtain 



R = (2cK5)§^§ 



(18) 



(19) 



and futhermore by use of R ■ T = const and eq. (17 b) 

1 



t = 



3 \ 2 ( To \ 2 



(20) 



2c 2 \koqoJ \T 

with the present day temperature Tq = 2, 7K of the background radiation. 



4.2 The condensation process 

In section 2 we realized that the cosmic matter gas runs from the end of recombination into 
a state very close to collision-dominated equilibrium. We describe this state in a very good 
approximation by the exact equilibrium as outlined in section 3. The single sharply edged 
gas clouds are then described by eqs. (15) where we consider the cosmic matter as a one- 
component gas with particle mass m p . Below we shall find that the temperature in this state 
is about 10 3 i^. Then we have even for r c ~ 10~ 5 (COBE) in the exponent of the density (15 
b) the enormous factor 



kTi 



(21) 
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Hence, when y drops only a very little below 1, say y = 0,999, the density falls practically 
to zero. 

Now we turn to the question how this process runs and how much time it takes. Clearly, the 
general relativistic Boltzmann equation 

*m-Ww =c{J) (22) 

determines this process. There is no hope to solve this equation unless we make some sim- 
plifications. Firstly we approcimate the collision term C(f) by the BGK-model [11-19 
] 

777 

C(/) = --^(/-/o) (23) 

where t c is the mean time between two collisions as in section 2 and /o is the distribution 
in collision-dominated equilibrium. Secondly we calculate the gravitational field T l - k only by 
means of the cosmological background (16). This means 

r^U=o = o 
rfa = l%s%, r? 4 = o (24) 
r o/3 = \ RR 5 al 3, r£ 4 = o = r| 4 

where a dot denotes derivative with respect to t and x = means x 1 = x 2 = x 3 = 0. We 
suppose that at the end of recombination small anisotropies of the background radiation 
are left from processes in the very early universe and consider the regions where the above 
normalized temperature variable y is below 1. In these regions the later (quasi)-equilibrium 
distribution /o approximately vanishes because the particle density is practically zero (see 
above). Futhermore we assume that in these regions 

f = f(t,p\p 2 ,p 3 ) (25) 

where the spatial momentum components p a are related with the 4-th component p 4 by the 
mass-shell condition 

9ijp¥ = g a( 3P a pP - (P 4 ) 2 = -m 2 p c 2 (26) 
In the cosmological epoch under consideration (T < 3 • 10 3 K, m = m p ) we have 

0a/?pV « (27) 

so that 

p 4 = m p c (28) 

Making use of the property (25) we can place the points we consider at the origin x = 0. 
Then the Boltzmann equation (22) reduces by means of the collision term (23) with /o = 
and the relations (24), (25) and (28) to 
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Next we integrate eq. (29) over the mass-shell (26) using the integration element [6] 

= \^~gd v x d v 2 dp 3 = \^d 3 P (30) 



where g = detgij. By means of the metric (16) we have 

7rU =0 = i? 3 (t)— (31) 
rripC 

By use of 

J f7r\x=o = n(t) = particle density (32) 

the integration yields 

9 f , r d f„3,^ d3 p 



m^ 



at J rripC J m p c 
|//vrU =0 -3|//vrU =0 



= n — 6—n 
R 



and 



f a 9f R 3 f a df 3 f 

J P W^ = ^cl P dp^ dP= ~ 3 J /7r '^° = 



(*) 

where we used in the integral (*) Gau' theorem. Hence we arrive at 

" = -(i + 4) n (33) 

The collision time t c is given by 

t c = (n ■ v t h ■ cry 1 (34a) 

with the termal velocity 

„„,= (—) (346) 

and the cross section 

a = irrl (34c) 
(r c = 7ao from eq. (2)). Futhermore it follows from eq. (19) that 

R 2 



R » (35) 
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Hence eq. (33) takes the form 



/'?--( AiT^n • - ) 



2 ( ?>k \ 2 

Ai = vrr c I 



(36) 



From eq. (20 we derive) 



i. 

T = T ( -i— ) 3 4 (37) 



Ac 4 K g J t a 

so that we obtain finally the differential equation 

/ A 2" 
n = — — n H — In 

Vt3 i, 



i 

A = 7rr 



\ m p I \4c 4 K Q 



(38) 



for the particle density n between the gas clouds. This density decreases on the one hand 
by collisions (first term in the paranthesis on the right hand side of eq. (38)) and on the 
other hand by the cosmological expansion (second term). Below we shall find that after 
recombination at first the collision term is greater than the expansion term by about eight 
orders of magnitude. However, this extreme superiority decreases very rapidly. 

The process (38) starts at the end of recombination at T = T re k = 3 • 10 3 K. Considering 
universes with 

go = a • 10 3 Vm~ 3 ; a = 0.5,1,2.5 (39) 
we obtain from eq. (20) the initial time 

t = t rek = • 0, 57 • 10 6 y (40) 
The particle density at t = t re k is given by n re k = Qrek/f^p where 

Qrek = a ■ Q rek , £ = 1, 4 • 10~ 21 £fCTO~ 3 (41) 

with g re k the density at the end of recombination in an universe with £o = 10~ 30 <7cm -3 . 
Hence we have the initial particle density 

n = n rek = a ■ 8, 2 • 10 2 cm" 3 (42) 

Next we write the differential equation (38) in a dimensionless form by introducing 

t = x- W 6 y, n = a ■ n ■ 8, 2 ■ 10 2 cm" 3 (43a) 

with the initial values 

Xrek = a~2 ■ 0, 57, ll rek = 1 (436) 



-9- 



Then eq. (38) reads 



11 = dx~ = + x )fi 



(44) 



f3 = al -8,0- 10 7 



The subsitution u = \i 1 leads to the linear differential equation 

2 



u = —u H r 

x X 3 



(45) 



with the general solution 



3 „ 2 



u{x) = Cx — -ftxz 



(46) 



The constant C is determined by the initial condition fi(x re k) = 1 Hence we arrive at 

.1-1 



H(x) = 



(47) 



For x well larger than x re k this solution can be written in the form 

(48) 



, x 4 1 



X 



J rek 



_ 4 
X 3 



The dimensionless form (44) of the differential equation (38) shows the enormous collision 
rate (first term in the paranthesies on the right hand side) just after recombination. 

We suppose that the process of cloud formation is worked out when the cosmological expan- 
sion rate predominaties more or less strongly the collision rate. Hence we determine the end 
x\ of the condensation process by 



-=b(3 

Xl 



1 ! 
"l 



6-10 



(49) 



which by means of the solution (48) yields 



Xl 



2, \^ 

-b + 1 J X rek 



From eq. (50) we obtain the cosmological instant 



(50) 



/j - ( h)+lY a ^ -0,57- 10 6 y (.11) 



of the end of the condensation process, its duration 

"2 



At 



6+1 -1 



0,57-a"2 -10 b y (52) 
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and by use of eq. (37) the temperature 



2 _1 
%rek \ 3 / 2 \ 2 



Finally we obtain from eq. (49) the value 



Aii = = 1 (54) 



6(16+1) 



2 A 



In tables 1 and 2 we give for b = 10 and 6 = 20 the values of ii, At, fj,i and Ti where a runs 
through the values according to eq. (39). 

The last row in table 1 gives the ratio ^cAt/R\ where R\ = R(t\) and R(ti) is calculated 
from eq. (19) by use of a mean Hubble-parameter h = 70. In section 4.3 we shall consider 
axially symmetric gas clouds. The maximal linear extension of these clouds is cAt because 
the clouds must be causally connected. From the centre of the clouds to the edge we have 
then a maximal extension ^cAt. This length must be much smaller than R± otherwise the 
weak field expansion (12) does not hold (see the line element (16)). 

4.3 The mass spectrum 

We can be sure that protogalaxies had been axially symmetric objects with their lower limit 
of spherically symmetric globular clusters. Hence we are interested in axially and spherically 
symmetric solutions of the single-cloud equations (15). 

We fix the particle mass as above by m = m p and let the anisotropies r c of the temperature 
vary in the range 

10~ 4 >t c > 10~ 7 (55) 

The background temperature T\ is determined by the parameter b (see tables 1 and 2). 
What about the density g c in the centre? In section 4.2 we established that the condensation 
process, that means the formation of the very sharp edges of the gas clouds, runs very rapidly. 
Furthermore the linear extension of the clouds is much smaller than the curvature radius R. 
Thus the clouds are not subject to the cosmic expansion and hence we can assume 

Qc = Qrek = UQ rek (56) 

with g re k and g re k according to eq. (41). As in the previous papersj [3,4] we introduce 
1 8vrG 



72 2 

Vefc C 



-Qrek, Irek = 6,6 ■ 10 5 ^ (57) 



and the dimensionless variable 



x a = —L- — x a (58) 

y Tclrek 



Then eqs. (15) take the form 
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A 5 y 



-oti exp 



(i-y) 



fcTi 

y<l, y(0,0,0) = l 



«i = -a, Q = aT 1 /g rek c 



(59a) 



£» = ae'refc exp 



fcTi 



where A 5 denotes the Laplacian with respect to x a . 



(596) 



4.3.1 Spherically symmetric clouds 

We denote the dinensionless radial coordinate by x so that y = y(x). Then eq. (59a) reads 

„2„ 



2 

y" + -y' = -ai exp 

x 



- g 



(60a) 



y(0) = 1 (606) 
To the natural boundary condition y(0) = 1 we add the regularity condition 

y'(0) = (60c) 

In his previous paper [3] the author investigated the differential equation (60) in detail. It 
turned out that the natural boundary condition (60b) implies the regularity (60c). The 
solution takes in the very vincinity of the centre the form 



ai + q a 1 +q . i 

H 3 — sin(pax) 

P^x 



P 



rripC r c «i 

m 



(61) 



This solution is used as start for a numerical integration of eq. (60a) by a Range-Kutta- 
method. 

The mass M of the gas clouds in the centre of the anisotropics is given by 



M 



= 47r / gr dr, r = ^JT c l rek x 
Jo 



(62) 



with the density g according to eq. (59 b). In section 4.2 we have shown that the particle 
density at the boundary of the clouds drops to ii\n rek when the condensation process is 
closed. Hence r\ in the mass integral (62) is strictly speaking given by g(r±) = [i\g c (see eq. 
(56)). However, the numerical calculations show that the masses practically do not change if 
we fix n by g(r±) = 10~ 5 £? c . This is due to the extremely rapid decrease of g. 

Writing the mass formula (62) in a dimensionless form we obtain 
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M 

M, 



- = 2, 13 • 10 18 • a ■ T C 2 / exp[- 7 (l - y(x))]x 2 dx 
Jo 



(63) 



7 



m p c 2 TC 
fcTi 



, exp[-7(l - = 10" 



The integral in eq. (63) can be easily evaluated by use of the numerical solution y(x). In 
tables 3 and 4 the results are listed where D = 2r± is the diameter of the clouds. All results 
are independent of the central value r c in the range (55). Finally we note that the radii \D 
of the clouds are very much smaller than R\ (see the end of section 4.2). 



4.3.2 Axially symmetric clouds 

We introduce dimensionless cylindrical coordinates x (distance from the symmetry axis) and 
C (z-coordinate) so that y = y(x,Q. Then eq. (59 a) reads 

d 2 y 1 dy d 2 y 

d^ + xd-x + dC = ~ ai eXP[ " 7(1 " V)] " q (640) 
with 7 according to eq. (63). Again we have the natural boundary condition 

2/(0,0) = 1 (646) 



and the regularity condition 



1(0.0 = o = |(.,o, 



(64c) 



Also the partial differential equation (64) has been investigated in detail by the author in his 
previous paper [3]. The solution is composed of an interior solution 



i{x,C) = 1 - j ( a i + 9) 



Cx 2 + (1 - C)C 



(65a) 



an exterior solution (when the density term on the right hand side of eq. (64 a) has practically 
fallen to zero) 



y a (x,() = 1 ~2 q 



\cx 2 +(i-c)e 



(656) 



and a fit function 



Uji.r.O- I \ < - v/2r(n -, + <,) \j hjx* + (1 - C)C 2 , 



e << 1 



(65c) 



with a flattening parameter C in the range 

< C < 1 



(65d) 



The fit function connects the interior and exterior solution. It is fit smoothly to the interior 
solution on a surface yi(x, Q = 1 — e = const, i.e. 



\cx 2 + (i - c)c 2 = -4- 

2 ai+q 



(66) 
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and intersects the exterior solution on the surface 

l C x 2 + {l-C)C 2 = %{^T+~q + V^i? (67) 
I q 

The solution is spherically symmetric for C = |. Then the interior solution (65 a) reads 

Vi(x,Q = I - ^ + q)(x 2 + C 2 ) (68) 

which agrees with the spherically symmetric solution (61) if one takes the first two terms in 
the expansion of the sin. 

On the surface of intersection (67) the density term on the right hand side of eq. (64 a) has 
the value 



«i exp[-7(l - y a )\ = ai exp 



— (V a i +Q + 



ei (69) 



We express the parameter e in terms of ei and obtain 

_ qln(ei/ai) 

7(V«i + 1 + \/«i) 2 



(70) 



Now we determine for given Ti(~> (7,7) and «i the parameter e by fixing e± such that the 
mass integral for C = | (spherical symmetry) agrees with the spherically symmetric masses 
calculated in section 4.3.1. As above we have the dimensionless mass formula 

XI Ci 

— = 2, 13 • 10 18 • a ■ t c 2 / / exp[- 7 (l - y(x, Q)]xdxdC (71a) 



where again at the upper limits of the integral the density has fallen to 10 -t, £ c ,ie. 





1-5, 



exp[- 7 (l-y(xi,0))] = 1(T 5 = exp[- 7 (l-y(x,Ci))], x < x x (716) 
and y(x, Q is the composed solution (65). 

In tables 5 and 6 the values of e\ are listed where we listed for reasons of later comparison in 
the last row the values of (i\ (see eq. (54)). Again the values of e\ do not depend on r c (the 
values of e, therefore, depend on r c since 7 depends on r c ). 

For C — > and C — ► 1 the clouds are pancake-like and cigar-like respectively (see fig. 2). We 
note that these configurations exist in an exact equilibrium state in the global rest frame (see 
section 3) , especially the pancakes do not rotate . They are hold by a complicated balance of 
pressure and gravitational forces where the gravitational forces come also from the radiation 
field [4]. 

What are the minimal and maximal value of C? We know from section 4.2 and eq. (56) that 
at the end of the condensation process the density at the boundary of the clouds drops to 
\i\Qc- Hence we obtain from eqs. (59 b), (56) and (58) a relation between the flattening 
parameter C and the semiaxises tq and zq by means of 
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exp[- 7 (l - y(x , 0))] = m = exp[- 7 (l - y(0, Co))] 



(72) 



x 



Co 



Tc^rek 



From tables 5 and 6 we see that e\/ai << fii. According to eq. (69) this means that the 
exponential is on the surface of intersection between the fit function and the exterior solution 
very much smaller than Hence y in eq. (72) lies in the domain of yf. By use of eqs. (65 
c) and (70) we obtain then from eq. (72) 



r VC 

trek 

piX -c) 



> = -0,96-l(T 5 (Ti[10 3 K])2 \ Inm + 



Irek 



qln(e 1 /a 1 ) 



V«i + q + ^/ai 



(V«i + q + v 7 "!) 2 J [-g(ai + g)in(ei/ai)]5 



(73) 



Since the clouds must be causally connected at the end of the condensation process the 
semiaxises ro for C — > and zq for C — > 1 are limited by 



2r ,max = cAt = 2z . 



(74) 



where At is the duration of the condensation process according to tables 1 and 2. With these 
upper limits of ro and zq we calculate from eq. (73) the values C m i n and C max . The results 
are listed in tables 7 and 8. 

With the composed solution (65) we can now easily compute from the mass formula (71) the 
masses of the clouds for given flattening parameters C where again it is sufficient to break 
off the integration when the density drops to 10 -5 £ c (instead of HiQ c )- The results are listed 
in tables 9 and 10. Again all is independent of the central value r c in the range (55). 

The mass spectrum we obtained is rather satisfactory. In fact, it shows the observed upper 
limit masses M ~ 10 12 M© as well as the lower limit masses M ~ 10 5 M Q where, beyond 
it, the lower limit masses are spheres. However, the upper limit masses M ~ 10 12 M Q are 
extremely thin pancakes, and galaxies are not observed in this form. Thus we investigate in 
the following section the evolution of the condensation " 'drops" ' where we argue as in the 
author's previous paper [3]. 

5) The evolution of the clouds after the condensation process 

After the condensation "drops" have formed a few 10 6 years after the big bang the universe 
cools further down, that means, the background temperature T\ decreases. If we assume that 

the quasi-equilibrium at first stands we conclude from eq. (73) that the semiaxises ro and zq 

i 

change by means of T x 2 in front of the curly bracket and q ~ Tf (see eq. (59a)). The second 
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term in the curly bracket can be neglected against ln[i\ and hence we have approximately 

1 1 



(75) 

l ± \ 



T 2 T 2 



This means that the density g d in the clouds decreases as 

Qd ~ T? (76) 
Then the characteristic expansion time of the clouds is given by 

,Qd, 2,Ti, 2R 



Tcxp W 9'Tx' 9R 



(77) 



where we used the cosmological law R T = const. We compare this characteristic time with 
the mean time t co u between two collisions of particles in the clouds. It is given by 

Tcoii = (n d -vth-a)' 1 (78a) 



n d = (786), 



v t h = (78c) 

\ m„ 



with, the particle density 
the thermal velocity 

and the cross section 

a = irr 2 c (78d) 
(see eq. (34)). According to eq. (76) the density q c \ can be written in the form 

Qd = Qd,i (t^) 1 ( 79fl ) 
where we use a mean value T\ = 10 3 K at the end of the condensation process and hence 

g dA ~ a ■ W~ 21 g ■ cm' 3 (79b) 
(see eq. (41)). In eq. (77) we use the field equation 

R 1/3X " (80) 



R c 2 \k q 



with the mean cosmological density g which is correlated to the background temperature T\ 
by means of g ~ Tf . Hence we have 

Wk) (81a) 

with the cosmological density g\ at T\ = 10 3 K given by 

Qi = (^) 3 a£refc = a ■ 5 • 1(T 23 gem- 3 (816) 
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Sticking all together we find 



^«-i.2.W^y (82) 
Tcxp \ T\ J 

This relation tells us that the clouds when expanding run rather rapidly out of their quasi- 
equilibrium state so that the contraction due to the decreasing temperature begins. We 
assume that the two effects nearly cancel so that the clouds after the condensation process 
at first simply drift apart until the universe has cooled down drastically, say down to T ~ 
100K = Ti- At this time the universe is about 10 8 y old and we assume that the clouds begin 
to contract. 

We consider this contraction in more detail. Especially we are interested in the upper limit 
masses M ~ 10 12 M Q which are extremely thin pancakes. Indeed, with C ~ 10 (see tables 
7 to 10) we obtain from eq. (73) 



C 1 (83) 



r Q y 2(1 - C) 4500 

How do these clouds contract? To answer this question we compare the gravitational energy 
Eq of the clouds with their thermal energy E t h after they have drifted apart and the universe 
has cooled down to T 2 ~ 100K. We have 

= B^P.V, P3 (84a) 

Va rn p 

with the volume V of the clouds and a nearly constant density 

g ~ oiQrek — ol ' 10~ 21 gcm~ 3 (84b) 
From these equations we obtain 

, , GrrugsMl/M^ , i/M\i 

k -WT 1 (d ~ 1 \~Mq) (85> 

which yields for the upper limit masses M ~ 1O 12 M 

\E G \ ~ 10 6 ^ (86) 

Hence we conclude that the semimajor axis of the upper limit clouds decrease in a rough 
approximation like a freely falling particle according to 

GM , . 

with the initial conditions 

r (0) = n , r(0) = (876) 

Here we use the initial value r(0) = r\ with ri according to eqs. (63) and (71) since the mass 
outside of T\ is negligible. The solution of eq. (87) reads 



V2GMt = r\ 
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With increasing contraction the cloud begins to rotate so that the contraction will be retarded 
until it stops at a final radius r 2 . We approximate the contraction time t con t r by 

tcontr = %tf (89) 

where tf is the free fall time (88) when r = \{r\ + r 2 ). 

We calculate the contraction times explicitly in an universe with qq = 1CT 30 gcm~ 3 (a = 1) 
for b = 10 (see table 9). For C = C m i n we have M = 6, 6 • 10 11 M Q which might be our galaxy 
including dark matter. We obtain the initial value r\ from eq. (73) by replacing [i\ by 1CP 5 
(see eq. (71)). We find 

n = 6, • W 5 ly (90) 
For the final value r 2 we take the present day value of our galaxy, hence 

r 2 = 0, 5 • 10 5 /y (91) 

Then we obtain a contraction time 

t contr = 2,5- 10 9 y (92) 

Calculating the red shift 

z = ^-l (93) 
corresponding to the time (92) from eq. (17) by use of a Hubble-parameter h = 70 we find 

z ~ 3 (94) 

Let us consider the same galaxy at the earlier time t = 0, 5t con t r - At this time the radial flow 
velocities v of matter are extremely high. From eq. (87) we obtain 



v(r) = V2GMJ- - — (95) 
y r n 

which for r = \{r\ + r 2 )(i = tf = 0, 5t contr ) yields 

v ~ 10 7 cm I sec (96) 
Comparing this velocity with the velocity of sound 



'kTY 2 



(97) 



for T ~ 300K we find 

v ~ 10 2 c s (98) 

Hence at the earlier time t = 0, 5t contr the galaxy is in a highly active state where the redshift 
corresponding to this time is 



z ~ 5,8 



(99) 
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Next we consider in table 9 the parameter C = 10 2 with a mass M = 10 M which is the 
upper limit of globular clusters. We have the ratio 

- = Ti ( 10 °) 
r 14 

and from eq. (85) 

\E G \ ~ 500^ (101) 

We calculate the contraction time as above where we assume that the contraction stops when 
the mass has become a sphere. Then we obtain 

t contr = l,6-10 8 y (102) 

However, doing the same for (7 = 0, 9999 with again M ~ 1O 7 M we find 

— = 71 (103) 

and a contraction time 

W = 3 • 10 9 y (104) 
Finally we obtain for C = 0, 999999 with a mass M ~ 10 s M the ratio 

— = 707 (105) 

and the enormous contraction time 

tcontr = 30 • 10 9 y (106) 

6) Conclusion 

Obviously our results agree very well with observations. The upper limit masses M ~ 1O 12 M 
evolve to fastly rotating spiral galaxies and are ready for the evolution of stars about 2, 5- 10 9 ?/ 
after the big bang. They are observed in their final state with a redshift z ~ 3 and in an 
earlier highly active state with a redshift z ~ 5. The lower limit masses M ~ 1O 5 M are 
born as spheres and are therefore ready for the evolution of stars a few 10 6 years after the 
big bang. Hence the stars in globular clusters are very old in agreement with observations. 

What about the cigar-like clouds for C — > 1? Maybe they have never existed. However, 
provided they have existed their contraction time is too long to evolve into clouds which are 
more or less spherically symmetric. They rather become unstable and decay into irregular 
clouds. Are these clouds the dark matter? 

Finally we discuss the fact that the masses of the clouds do not depend on the central 
value t c in the range (55) (see section 4.3). First of all we note that the anisotropies must be 
axially symmetric only in the very inner region since the density of the clouds falls practically 
to zero when the variable y = t/t c drops to y = 0, 99 (see eq. (59 b)). Hence an anisotropy 
r c ~ 10~ 5 can be composed of say 100 axially symmetric anisotropies r c ~ 10~ 7 and covers 
then about 100 galaxies. Thus the spatial pattern of anisotropies of the background radiation 
after recombination might have been composed of chains and slices and hence matter has been 
arranged in chains and slices of galaxies. To this initial pattern of galaxies the machinery 
used by the Virgo Consortium [2] can be applied to give the large scale structures we observe 
today. However, in contrast to the results by the Virgo Consortium the structures consist 
now of luminous matter. 



-19- 



REFERENCES 

[I] S. Weinberg, Gravitation and Cosmology, 
J. Wiley, New York (1972) 

[2] Virgo Consortium (F. Pearce, P. Thomas, A. Jenkins, 
C. S. Frenk, H. Couchman, S. White, J. Colberg, 
G. Efstathiou, A. Nelson, J. Peacock), 
Annual Report 1995, Max-Planck-Institute Garching 

[3] G. Lessner, GRG 26 (1994), 385 

[4] G. Lessner, GRG 27 (1995), 417 

[5] G. Lessner, Astrophys. and Space Sci. 161 (1989), 175 

[6] J. M. Stuart, Cargse Lectures in Physics 6, 

E. Schatzman, ed., Gordon and Breach, New York (1973) 

[7] W. Israel, in General Relativity, L. O'Raifeartaigh, ed. 
Clarendon Press, Oxford (1972) 

[8] H. Dehnen, O. Obregon, Astron. and Astrophys. 12 (1971), 161 

[9] B. Chaboyer, Y.C. Kim, Ap. J. 454 (1995), 767 

[10] M. Salaris, S. DegPInnocenti, A. Weiss, Ap. J. 479 (1997), 665 

[II] S. Chapman, T.G. Cowling, The Mathematical Theory of 
Non-Uniform Gases, Cambridge University Press, Cambridge (1970) 

[12] C. Cercignani, Theory and Application of the Boltzmann Equation, 
Scottish Academic Press, Edinburgh, London (1975) 

[13] C.Marie, Ann. Inst. H. Poincare 10 (1969), 67 

[14] J. D. Nightingale, Ap. J. 185 (1973), 105 

[15] J. L. Anderson, H. R. Witting, Physica 74 (1974), 466 

[16] J. L. Anderson, A. C. Payne, Physica 85 A (1976), 261 

[17] R. Dominguez-Tenreiro, R. Hakim, Phys. Rev. D 15 (1977), 1435 

[18] A. Majorana, J. Math. Phys. 31 (1990), 2042 

[19] F. P. Wolvaardt, R. Maartens, Class. Quant. Grav. 11 (1994), 203; 
14(1997), 535 



-20- 



a 


0.5 


1 


2 


5 


*i[10 6 y] 


3,7 


2,6 


1,9 


1,2 


At[10 6 ?/] 


2,9 


2,1 


1,5 


0,92 


lii 


1,8 -KT 9 


1,3 -KT 9 


9,2 -10- iU 


5,8 -10- iU 


^cAt/Ri 


0,04 


0,03 


0,02 


0,01 



Table 1: The values of ti, At and /ii 

for b = 10 (Ti = 1,08- 10 3 iT) 



a 


0.5 


1 


2 


5 


*i[10 e y] 


5,9 


4,2 


3,0 


1,9 


At[10 6 ?/] 


5,1 


3,6 


2,6 


1,6 


/ii 


6,7 -10- 1U 


4,8 -10- 1U 


3,4- 10" 1U 


2,1 -10- 1U 



Table 2: The values of ti, At and /ii 

forb = 20 (Ti = 0,792- 10 3 K) 
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a 


0.5 


1 


2 


5 


M/Mq 


2,7 -10 5 


2,5 -10 5 


2,4 -10 5 


2,2 -10 5 


D[iy] 


870 


820 


770 


700 



Table 3: Masses M and diameters D of gas 
clouds in the centre if spherically 
symmetric anisotropics for b = 10 
(Ti = 1, 08 -10 3 i^ = 0,0083) 



a 


0.5 


1 


2 


5 


M/Mq 


2,8 -10 5 


2,6 -10 5 


2,5 -10 5 


2,4 -10 5 


D[iy] 


1240 


1160 


1070 


960 



Table 4: Masses M and diameters D of gas 
clouds in the centre of sherically 
symmetric anisotropics for b = 20 
(7i = 0, 792 • 10 3 if, q = 0, 0024) 
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a 


0.5 


1 


2 


5 


£1 


6 -lO" 16 


2 -10" 24 


1Q -38 


10" 7i 




1,8 TO" 9 


1,3 -10" 9 


9,2 -10- 1U 


5,8 -10- 1U 



Table 5: The parameter s\ for b = 10 
(Ti = 1,08- 10 3 if, q = 0,0083) 



a 


0.5 


1 


2 


5 


£i 


1Q -35 


1Q -56 


5 • 10" 9U 


10 -164 


111 


6,7 -10- 1U 


4,8 -10- 1U 


3,4 -10- 1U 


2,1 -10- 1U 



Table 6: The parameter e\ for b = 20 

(Ti = 0, 792 • 10 3 if, q = 0,0024) 
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a 


0.5 


1 


2 


5 


C ■ 


1,21 • 10" 7 


1,50- 10" 7 


1,87- 10" 7 


2,75- 10" 7 


Cmax 


0,999999 
939 


0,999999 
925 


0,999999 
906 


0,999999 
862 


Table 7: The values of C m ^ n 
and Cmax for b = 10 
(Ti = 1,08- 10 3 K, q = 


= 0,0083) 




a 


0.5 


1 


2 


5 


C ■ 


4,66- 10" 8 


5,95- 10" 8 


7, 34 • 10" 8 


1,10- 10" 7 


Cmax 


0,999999 
976 


0,999999 
970 


0,999999 
963 


0,999999 
945 



Table 8: The values of C m j n 
and Cmax for b = 20 
(Ti = 0, 792 • 10 3 K, q = 0,0024) 



O 


M / Mq^Q! — U, 0) 


A/f / A/f 
iVi / IVIq 


[a - 1) 


A/f / A/f ^/ o\ 

M / 1V1q[CX — Z) 


A/f / A/f 
IVl J IVIq 


( ~. r:\ 

\0L — o) 




or: i nil 

o,o •1U 


p p 
0,0 


1 nil 
1U 


a n 1 nil 

4,y -iu 


9 1 

3,1 • 


1 nil 
IU 


1 n— 6 
1U 


1 n 1 nil 
l.U -IU 


n o 

y,o 


1 nlO 
1U 


n 1 1 nlO 

y,i -IU 


o,o • 


1 nlO 
IU 


1 n — 5 
1U 


1 n 1 nlO 
l.U -IU 


n o 

y,o 


•1U 


n 1 1 n9 

y,i -IU 


o,0 


1 n9 
•IU 


1 n— 4 
1U 


1 n 1 n9 
l.U -IU 


n o 

y,o 


•1U 


n 1 1 n8 

y,i -IU 


o,0 


1 n8 
•IU 


1 C\— 3 
1U 


1 n 1 n8 
l.U -IU 


n o 

y,o 


1 n7 
•1U 


n 1 1 n7 

y,i -IU 


0,0 


1 n7 
•IU 


1 n-2 
1U 


1 n 1 n7 
l.U -IU 


n o 

y,o 


1 n6 
•1U 


no 1 n6 

y,z -iu 


O P 

0,0 


1 n6 
•IU 


n 1 
U, 1 


1 n 1 n6 
l.U -IU 


1 n 
1,0 


1 n6 
•1U 


n £ 1 n5 

y,o -IU 


n n 

y,u 


1 n5 
•IU 


1 /Q 

z/o 


o i 1 n5 
z, / -IU 


z,0 


1 n5 
•1U 


o /i 1 n5 
z,4 -IU 


o o 

z,z 


1 n5 
•IU 


u, y 


o /? i n5 
O,0 -1U 


Q /I 

3,4 


1 n5 
•1U 


o o 1 n5 
o,z -IU 


q n 
o,U 


1 n5 
•IU 




1 n -in 6 


q q 


•1 n 5 


q 9 -1 n 5 

<J,Z; _LU 


8 6 


•1 n 5 


0, 9999 


1.0 -10 7 


9,8 


•10 6 


9,1 -10 6 


8,5 


•10 6 


0, 999999 


1.0 -10 8 


9,8 


•10 7 


9,1 -10 7 


8,5 


•10 7 


Cmax 


4.2 -10 8 


3,6 


•10 8 


3,0 -10 8 


2,3 


•10 8 



Table 9: Masses of gas clouds in the centre 
of axially symmetric anisotropics 
for b = 10 (Ti = 1,08 • 10 3 K, q = 



0,0083) 



U 


l\/T 1 l\/T ( y^, n £\ 
M / IVIq [a = U, D) 


M / 1VLq\ol =1; 


1\/T 1 1\/T 

IVl / IVIq 


( a - 2 ) 


71 /f / ft /f 

M /IVIq 


(a = oj 




z,d -1U 


1 7 1 nl2 
1,/ -1U 


1 9 

1,3 • 


1 nl2 
1U 


Q P. 

o,o • 


1 nil 
1U 


1 n-7 
1U 


1 1 ml2 
1.1 -1U 


1 n 1 r>12 
1,U •11) 


n i 

y,/ • 


1 nil 
1U 






1 n— 6 
1U 


11 1 nil 

1.1 -1U 


1 n 1 nil 
1,U -1U 


n i 

y, / • 


1 nlO 
1U 


n /i 

y,4 • 


1 nlO 
1U 


1 n— 5 
1U 


11 1 nlO 
1.1 -1U 


1 n 1 nlO 
1,U -1U 


n i 

9,7 


1 n9 
•1U 


n /i 

y,4 


1 n9 
•1U 


1 n— 4 
1U 


11 1 n9 
1.1 -1U 


1 n 1 n9 
1,U -1U 


n 7 

9,7 


1 f|8 
•1U 


n /i 

y,4 


1 n8 
•1U 


1U 


11 1 n8 
1.1 -1U 


1 n 1 n8 
1,U -1U 


9,7 


1 n7 
•1U 


n /i 

y,4 


1 n7 
•1U 


1 n-2 
1U 


11 1 n7 
1.1 -1U 


1 n 1 n7 
1,U -1U 


n o 

y,o 


1 n6 
•1U 


n c 

y,o 


1 n6 
•1U 


U, 1 


11 1 n6 
1.1 -1U 


11 1 n6 
1,1 -1U 


1 n 
1,0 


1 n6 
•1U 


y,y 


1 n5 
•1U 


Z/6 


o o 1 n5 
z,o -1U 


o a 1 n5 
Z,0 -1U 


Z,0 


1 n5 
•1U 


Z,0 


1 n5 
•1U 


n o 

u, y 


o,o -1U 


o,b -1U 


9 /I 

3,4 


1 n5 
•1U 


9 9 

o,o 


1 n5 
•1U 




1 1 -1D 6 


1 -10 6 


Q 8 


•1 n 5 


Q 5 


• 10 5 


0, 9999 


1.1 -10 7 


1,0 -10 7 


9,7 


•10 6 


8,5 


•10 6 


0, 999999 


1.0 40 8 


1,0 -10 8 


9,7 


•10 7 


9,4 


•10 6 


Cmax 


7.0 -10 8 


5,9 -10 8 


5,1 


•10 8 


4,0 


•10 8 



Table 10: Masses of gas clouds in the centre 
of axially symmetric anisotropics 
for b = 20 (Ti = 0, 792 • 10 3 K, q 



= 0,0024) 



tu [10 V] 




-30 -3 n 

P o [l0 gem ] 



1 2 3 4 5 6 7 8 9 10 11 12 13 



Fig. 1: The age t u of the universe depending 
on the Hubble-parameter 
H D = h -km/sec ■ Mpc and the 
present day density P D 



